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(57) Abstract 



A method and apparatus for controlling a physical system (10) responsive to an input to produce a selected output comprises defining 
a virtual identity system that includes the physical system for receiving the input to provide an actual output. A quality of identity of the 
virtual identity system is checked using at least a function of the actual output. In a preferred embodiment, adjustments are made to a 
model of the physical system as a function of the quality of identity in order to render the correct input more effciently and accurately. 
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METHOD AND APPARATUS FOR GENERATING INPUT 

SIGNALS IN A PHYSICAL SYSTEM 

BACKGROUND OF THE INVENTION 
The present invention relates to a control 
5 of a system, machine or process that is repetitive in 
nature or is amenable to at least some degree of 
rehearsal. More particularly, the present invention 
relates to calculating a model to be used for 
generating drive signals as input to a vibration 
10 system.. 

Vibration systems that are capable of 
simulating loads and/or motions applied to test 
specimens are generally known. Vibration systems are 
widely used for performance evaluation, durability 
15 tests, and various other purposes as they are highly 
effective in the development of products. For 
instance, it is quite common in the development" of 
automobiles^ motorcycles, or the like, to subject the 
vehicle or a substructure .thereof to a laboratory 
20^- - environment- thaf simulates operating conditions such 
as a road or test track. Physical simulation • in the 
laboratory involves a well-known method of data 
acquisition and analysis in order to develop drive 
signals that can be applied to the vibration system to 
25 reproduce the operating environment. This method 
includes instrumenting the vehicle with transducers 
"remote" to the physical inputs of the operating 
environment. Common remote transducers include, but 
are not limited to, strain gauges, accelerometers , and 
30 displacement sensors, which implicitly define the 
operating environment of interest. The vehicle is then 
driven in the same operating environment, while remote 
transducer responses (internal loads and/or motions) 
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are recorded. During simulation with the vehicle 
mounted to the vibration system, actuators of the 
-vibration system ar-e driven so as to reproduce the 
recorded remote transducer responses on the vehicle in 
5 the laboratory. 

However, before simulated testing can. occur, 
the relationship between the input drive signals to 
the vibration system and the responses of the remote 
transducers must be characterized in the laboratory. 
10 Typically, this "system identification" procedure 
involves obtaining a respective model or transfer 
function of the complete physical system (e.g. 
vibration system, test' specimen, and remote 
transducers) hereinafter referred to as the "physical 
15 system" ; calculating an inverse model or transfer 
function of the same; and using the inverse model or 
transfer function to iteratively obtain suitable drive 
signals for the vibration system to obtain 
substantial! y^-- the skme response from the remote 
20 transducers on the test specimen in the laboratory 
situation as v/as found in the operating environment. 

As those -skilled in the art would 
appreciate, this process of obtaining suitable d.rive 
signals is not altered when the remote transducers are 
25 not physically remote from the test system inputs 
(e.g. the case where "remote" transducers , are the 
feedback variables, such as force or motion, of the 
■vibration system controller) . 

Although the above -described system and 
3 0 method for obtaining drive signals for a vibration 
system has enjoyed substantial success, there is a 
continuing need to improve such systems . In 
particular, there is a need to improve models of the 



SUBSTITUTE SHEET (RULE 26) 



( 



wo 99/38054 



r 



PCT/US99/01233 



physical system and the iterative process for 
obtaining the drive signals. 

SUMMARY OF THE IIsTVENTION 
An aspect of the present invention relates to a 
5 method and apparatus or . system controller for 
controlling a physical system responsive to an input 
to produce a selected output . The method comprises 
defining a virtual identity system that includes the 
physical system for receiving the input to provide an 
10 actual output; and checking a quality of identity of 
the virtual identity system using at least a function 
of the actual output. The system controller includes 
program modules to perform the method. Instructions 
can also be provided on a computer readable medium, to 
15 implement the method. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of an exemplary 
environment for practicing the present invention. 

FIG. 2 is a computer for implementing the 
2 0 present invention. 

FIG. 3A is a flow chart illustrating the 
steps involved in an identification phase of a prior 
art method of vibration testing. 

FIG. 3B is a flow chart illustrating the 
25 steps involved in. an iterative phase of a prior art' 
method of vibration testing. 

FIG. 3C is a flow chart illustrating the 
steps involved in another iterative phase of a .prior 
art method of vibration testing. 
30 FIG. 4A is a detailed block diagram of a 

prior art . iterative process ror obtaining drive 
signals for a vibration system with an adjuster of the 
present invention. 
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FIG. 4B is a detailed block diagram of 
another prior art iterative process for obtaining 
-drive signals for a -vibration system with the adjuster 
of the present invention. 
5 FIG. 5 is a general block diagram of an 

aspect of the present invention. 

FIG. 6 is a detailed block diagram of an 
embodiment of the invention of FIG. 5. 

FIG. 7 is a flow chart illustrating the 

10 steps involved for operating the embodiment of FIG. 6. 

FIG. 8 is a detailed block diagram of 
another embodiment of the invention of FIG. 5. 

FIG. 9 is a general block diagram of a 
second aspect of the present invention. 
15 FIG. 10 is a general block diagram of a 

third aspect of the present invention. 

FIG. 11 is a detailed block diagram of an 
embodiment of the invention of FIG. 9. 

FIG. 12 is a detailed block diagram c-^ an 
20 embodiment of the invention of FIG. 10. 

FIG. 13 is a block diagram of a fourth 
aspect of the present invention. 

FIG. 14 is a block diagram. of a fifth aspect 
of the present invention. 
25 FIG. 15 is a block diagram of a sixth aspect 

of the present invention. 

FIG. 16 is a: pictorial representation of a 
seventh aspect of the present invention. 

FIG. 17 is a block diagram of a eighth 
30 aspect of the present invention. 

FIG. 18 is a block diagram of ninth aspect "of the 
present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

FIG. 1 illustrates a physical system 10. The 
physical system 10- generally includes a vibration 
system 13 comprising a ser"o controller 14 and an 
5 actuator 15. In the schematic illustration of FIG. 1, 
the actuator 15 represents one or more actuators that 
are coupled through a suitable mechanical interface 16 
to a test specimen 18 . The servo controller 14 
•provides an actuator command signal 19 to the actuator 
10 15, which in turn, excites the test specimen 18. 
Suitable feedback 15A is provided from the actuator 15 
to the servo controller 14 . One or more remote 
transducers 20 on the test specimen 18, such as 
displacement sensors, strain gauges, accelerometers , 
15 or the like, provide a measured or actual response 21. 
A physical system controller 23 receives the actual 
response 21 as feedback to compute a drive 17 as input 
to the physical system 10, In an iterative process 
d-is-vjussed below, the physical system controller 23 
2 0 generates the drive 17 for the physical system 10 
based on the comparison of a desired response provided 
at 22 and the actual response 21 of the remote 
transducer 20 on the test specimen. 18. Although 
illustrated in FIG. 1 for the single channel case, 
25 multiple channel embodiments with response 21 
comprising N response components and the drive 17 
comprising M drive components are typical and 
considered another embodiment of the present 
invention . 

30 Although described herein where the physical 

system comprises the vibration system 13 and remote 
transducer 20, aspects of the present invention 
described below can be applied to . other physical 
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systems . For instance, in a manufacturing process, the 
physical system includes the manufacturing machines 
(-e.g. presses, molding apparatus, forming machines, 
etc.) and the drive 17 provides command signals to 
5 said machines, and the actual response 21 comprises 
manual or automatic measured parameters of the 
manufactured article- such as a critical dimension. 
Another example includes an oil refinery where the 
physical system is the process plant and the actual 

10 response 21 comprises intermediate or final parameters 
related to output products. 

FIG. 2 and the related discussion provide a 
brief, general description of a suitable computing 
environment in which the invention may be implemented . 

15 Although not required, the physical system controller 
23 will be described, at least in part, in the general 
context of computer-executable instructions, such as 
program modules, being executed by a computer 30. 
Generally, program modules include routine programs, 

2 0 objects, _components^, data structures,^ etc., which 
perform particular tasks or implement particular 
abstract data types. The program modules are 
illustrated below using block diagrams and flowcharts. 
Those skilled in the art can implement the block 

2 5 diagrams and flowcharts to computer-executabJLe 

instructions. Moreover, those skilled in the art will 
appreciate that the invention may be practiced with 
other computer system configurations, including multi- 
processor systems, networked personal computers, mini 

3 0 computers, main frame computers^ and the like. The 

invention- may also be practiced in distributed 
computing environments where tasks are performed by 
remote processing devices that are linked through a 
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communi cat ions network. In a distributed computer 
environment, program modules may be located in both 
-local and remote memory storage devices . 

The computer 30 illustrated in FIG. 2 
5 comprises a conventional personal ■ or desktop computer 
having a central processing unit (CPU) 32, memory 34 
and a system bus 36, which couples various system 
components, including the memory 34 to the CPU 32. The 
system bus 3 6 may be any of several types of bus 
10 structures including a memory bus or a memory 
controller, a peripheral bus, and a local bus using 
any of a variety of bus architectures. The memory 34 
includes read only memory (ROM) and random access 
memory (RAM) , A basic input/output (BIOS) containing 
15 the basic routine that helps to transfer information 
between elements within the computer 30, such as 
during start-up, is stored in ROM. Storage devices 38, 
such as a hard disk, a floppy disk drive, an optical 
disk drive, etc., are coupled to the system bus 36 and 
20 are used for storage of ■ programs and data. It should 
be appreciated by those skilled in the art that other 
types of computer readable media that are accessible 
by a computer, such as magnetic cassettes, flash 
memory cards, digital video disks, random access 
25 memories, read only memories, and the like, may also 
be used as storage devices. Commonly, programs are 
loaded into memory 34 from at least one of the storage 
devices 3 8 with or without accompanying data. 

An input device 4 0 such as a keyboard, 
0 pointing device (mouse) , or the like, allows the user 
to provide commands to the computer 30. A monitor 42 
or other type of output device is further connected to 
the system bus 36 via a suitable interface and 
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provides feedback to the user. The desired response 22 
can be provided as an input: to the computer 3 0 through 
- H communications li-nk, such as a modem, or through the 
removable media of the storage devices 38. The drive 
5 signals 17 are provided to the physical system 10 of 
FIG. 1 based on program m.odules executed by the 
computer 30 and chrough a suitable interface 44 
coupling the computer 30 to the vibration system 13. 
The interface 44 also receives the actual response 21. 

Before describing the present invention, it 
may also be helpful no review, in detail, a known 
method for modeling che physical system 10 and 
obtaining the drive 17 to be applied thereto. Although 
described below with respect to a test vehicle,, it 
15 should be understood that this prior art m.ethod and 
the present invention discussed belov/ are not confined 
to testing only vehicles, but can be used on other 
types of test specimens and substructures or 
cor.poiients thereof. In addition, the descriptio.n is 
20 <^onB assuming spectral analysis ^ based modeling 
estimation and implementation though operations can be 
carried by several other mathematical techniques (e.g. 
Adaptive Inverse' Control (AIC) type models-, parametric 
regression techniques such as Auto Regressive 
25 Exogenous (ARX) and . State Space types of models, or 
combinations thereof) . 

Referring to FIG. 3A, at step 52, the cest 
vehicle is instrumented • v^ith the remote transducers 
20. At step -54, the vehicle is subjected to the field 
3 0 operating environment of interest and the remote 
transducer responses are measured and recorded. For 
instance, the vehicle can be driven on a road or test 
track. The measured remote transducer responses. 
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typically analog, are stored in the computer 3 0 in a 
digital format through analog- to-digital converters, 
-as is commonly known. 

Next, in an identification phase, the 
5 input/output model of the physical system 10 is 
determined. This procedure includes providing drive 17 
as an input to the physical system 10 and measuring 
the remote transducer response 21 as an output at step 
56. The drive 17 used for model estimation can be 
10 random "white noise" having frequency components over 
a selected bandv/idth. At step 58, an estimate of the 
model of the physical system 10 is calculated based on 
the input drive applied and the remote transducer 
response obtained at step 56. In one embodiment, this 
15 is commonly known as the "frequency response function" 
(FRF) . Mathematically, the FRF is a N x M matrix 
wherein each element is a frequency dependent complex 
variable (gain and phase versus frequency) . The 
columns of the matrix correspond to the inputs, while 
20 the rows correspond to the outputs. As app_reciated by 
those skilled in the art, the FRF may also be obtained 
directly from prior tests using the physical system 10 
or other systems substantially similar to the physical 
system 10 

25 An inverse model H(f)'^ is needed to 

determine the physical drive 17 as a function of the 
remote responses at step 60. As appreciated by those 
skilled in the art, the inverse model can be 
calculated directly. Also, the , term "inverse" model as 

30 used herein includes a M x N "pseudo- inverse " model 
for a non- square N x M system. 

At this point in the prior art, the method 
enters an iterative phase, illustrated in FIGS. 3B and 
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4A, to obtain drive 17 which produces actual response 
21 that ideally replicates the desired remote 
-transducer response 22 (hereinafter "desired 
response"). The inverse physical system model H(f)'^ is 
5 represented at 12, while physical system (vibration 
system, test vehicle, remote transducers and 
instrumentation) is represented at 10. Referring to 
FIG. 3B, at step 78, the inverse model 72. is applied 
to a target response correction 77 in order to 
LO determine an initial drive 17 Xi{t). The target 
response correction 77 can be the desired response 22 
for the initial drive, though most often it is reduced 
by a relaxation gain factor 95. The calculated drive 
17 Xi(t) from the inverse model 72 is then applied to 
.5 the physical- system 10 at step 80. The actual remote 
transducer response 21 (hereinafter "actual response") 
y^Ct) of the physical system 10 to the applied drive 17 
x^(t) is then obtained at step 86. If the complete 
physical system 10 is linear (allowing a relaxation 
0 gain ^95- of unity) , then the initial drive 17 Xi(t) 
could be used as the required drive. However, since 
■physical systems are typically non-linear, the correct 
drive 17 has to be arrived at by an iterative process. 
(As appreciated by those skilled in the art, drive 17 
5 used in previous tests for a similar physical system 
may be used as the initial drive.) 

The iterative process involves recording ^the 
first actual response yi(t) resulting from the initial 
drive Xi(t) and comparing it with the desired response 
0 22 and calculating a response error 89 Ay^^ as the 
difference at step 88. (The first actual response 
signal yi(t) is provided at 87 in FIG. 4A. ) The 
response error 8 9 Ay^ is compared to a preselected 



SUBSTITUTE SHEET (RULE 26) 



r 



wo 99/38054 PCT/US99/01 233 

-11- 

threshold at step 90 and if the response error 89 
exceeds the threshold an iteration is performed. 
"Specifically the response error 89 Ay^ is reduced by 
the relaxation gain factor 95 to provide the new 
5 target- response correction 77. In this embodiment, the 
inverse transfer function H(f)*^ is applied to the new 
target .response correction 77 to create a drive 
correction Axj 94 (step 91) that is added to the first 
drive y:^{t) 17A to give a second drive X2(t) 17 at step 
10 92. The iteration process (steps 80-92) is repeated 
until the response error 89 is brought down below the 
preselected threshold on all channels of the response. 
The last drive 17, which- produced a response 21, that 
was within the predetermined threshold of the desired 
15 response 22, can then be used to perform specimen 
■ testing. 

As described, the response error 89 Ay is 
commonly reduced by the relaxation gain factor (or 
iteration gain) 95 .to form the target response 

'20 corf ect iorT 7 7T~ The iteratxon gaTn" '9"5 ' stabilizes the 
iterative process and trades off rate-of -convergence 
against iteration overshoot. Furthermore, the 
iteration gain 95 minimizes the possibility that the 
test vehicle will be overloaded during the iteration 

25 process due to non-linearities present in the physical 
system 10. As appreciated by those skilled in the art, 
an iteration gain can be applied to the drive 
correction 94 Ax and/or the response error 89. It 
should be noted in FIG. 4 A that storage devices 3 8 can 

30 be used to store the desired response 22, the actual 
responses 21 and previous drives 17A during the 
iterative process. Of course, memory 34 can also be 
used. Also, a dashed line 93 indicates that the 
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inverse model 72 is an estimate of the inverse of the 
physical system 10. The block diagram of FIG. 4A, as 
-discussed above, can be implemented by those skilled 
in the art using cc"nmercially available software 
5 modules such as included with RPCIII from MTS Systems 

Corporation of Eden Prairie, Minnesota. 

At this point, a modified method of the 
- prior ■ art for calculating the drive can also be 
discussed. ■ The modified prior' art method 'includes the 
10 steps of the identification phase illustrated in FIG. 
3A and many of the steps of the iterative phase 
illustrated in FIG. 3B. For convenience, the iterative 
steps of the modified method are illustrated in FIG. 
3C and the block diagram as illustrated in FIG. 4B: As 
15 illustrated in FIG. 4B, the calculation of the target 
response correction 77 is identical. However, if the 
response error 89 between the actual response 21 and 
the desired response 22 is greater than a selected 
threshold, then the target response correction 77 is 

20- ...added -to a... previous . target, response 79A-at-.step 97-. to . 

obtain a new target response 79 for the current 
iteration. The inverse model 72 is applied' to the 
target response 79 to obcain the new drive 17. As 
illustrated in FIG. 4B, the iteration gain 95 can be 
25 .used for the reasons discussed above. 

Generally, an aspect of the . present 
invention includes the response time history error 
iteration loop described above with respect to FIGS. 
4A and 4B, while including an adjuster 100 that 
30 operates during each step of the iterative process, to 
improve the physical system inverse model 72 . As 
illustrated in FIG. 4A, the adjuster 100 corrects the 
inverse model 72 which receives the' target response 
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correction 77 directly as a simple function of the 
response error 8 9 (i.e. without previous target 
-information 79A of - FIG. -4B) and where the physical 
system drive 17 comprises drive' correction 94 in 
5 combination with a previous drive 17A. Conversely, as 
illustrated in FIG. 4B, the inverse model 72 receives 
the target response 7 9 as the combination of the 
target response correction 77 and a the previous 
.^target' response- 79A., and drive 17 is directly obtained 
10 by applying the inverse model 72. In the case of . FIG 
4B, the adjuster 100 corrects the inverse model 72 in 
a conceptually identical fashion as in FIG. 4A'. 
However, as will be discussed below, the 
configurations of FIGS. 4A and 4B render different 
15 signals available to the virtual identity modeling 
process each with inherent situational advantages. 
Furthermore, as will be described below, the adjuster 
100 can also operate in an iterative manner. 

An aspect of the present invention is 

_ 2 0 illustrated in the block _cy.agram of _ FIG 5 Generally, 

this aspect of the invention includes a method of 
controlling the physical system - 10 to produce an 
actual response 21 that ideally matches the desired 
response 22 consistent with prior art as discussed. 
2 5 The method includes generating an inverse model _98 
(e.g. the inverse transfer function H(f)'^) of - the 
physical system 10, wherein the inverse model 9 8 is 
applied to the target response 79 to obtain the drive 
17 intended to generate the desired response 22 from 
30 the physical system 10 as the actual response 21. As 
will be described below with respect to FIG. 6, the 
target response correction 77 can be combined with the 
previous target response 79A, when the inverse model 
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98 is applied to the target response 79 to realize the 
complete drive 17. Alternatively, as illustrated in 
-FIG. 8 with dashed -lines/ the inverse model 98 can be 
applied to the target response correction ' 77 directly 
5 to realize the drive correction 94, which is then 
subsequently combined with the previous drive 17A to 
provide a new drive 17 for the physical system 10. 

Referring back to FIG. 5, the adjuster 100 
can comprise a virtual identity system, modeler that 

10 checks a quality of identity of the inverse model 98 
in combination with the physical system 10. The 
combination of the inverse model 98 and the physical 
system 10 is designated as a virtual identity system 
96. The quality of identity, as measured by .the 

15 virtual identity system modeler 100, assesses the 
accuracy of the physical system 10 model relative to 
the operating characteristics of the physical system 
10. In this embodiment, the quality of identity is 
measured via the series connection of the physical 

20 system 10 with_ the_ inverse system model 98 . When the 
inverse system model 98 identically matches the 
physical system 10, the quality of identity 
calculation produces an identity result, indicating an 
ideal inverse system model estimate. 

2 5 The advantage of this approach is that , a 

model of a potential correction to the physical system 
inverse model that improves the quality of the virtual 
identity system can take a simpler form than the 
inverse model itself. Consequently, the model 

3 0 correction is much easier to obtain than a new 

estimate of the complete inverse system model. This 
approach leads to a straightforward algorithmic 
formulation of such a correction model estimate (e.g. 
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the FRF between the target response and the actual 
response) . The simple form of this modeling technique 
a.llows smaller segments of noisier and more correlated 
data to be used, thereby providing an advantage in an 

5 adaptive environment for either non -parametric 
(spectrum analysis, etc.) or parametric (ARX, etc.) 
modeling methods. 

In one embodiment of FIG. 5, with switches 
lOlA and lOlB at position "G" , the virtual identity 

0 system modeler 100 compares the target response 7 9 and 
the actual response 21. In another embodiment with 
switches lOlA and lOlB at position "L" , the virtual 
identity system modeler 100 compares the target 
response correction 77 and an actual response 

5 correction 103 (difference between response y^^ and 

response yi^i) obtained from application of the drives 

17 to the physical system 10 for successive 
iterations. In either embodiment, the virtual identity 
• system modeler 10 0 subsequently adapts the inverse 
model 98 as a function of the quality of identity from 
iteration to iteration. In other words, the virtual 
identity system modeler 100 provides as an output, 
model (e.g. VrF) correction values 104 to adjust at 
least some of the values present in the inverse model 
98. In both of these preferred embodiments, the 
correction values 104 are derived on a channel by 
channel basis, a natural simplification facilitated by 
the virtual identity system. Cross -coupling effects do 
not then form part of the correction values 104, 
although all terms of the inverse model 98 will 
generally change when the correction values 104 are 
applied. These preferred embodiments however do not 
preclude' the option of including some or all cross- 
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coupling terms into the correction values 104. 

It should be understood that switches lOlA 
-and lOlB represent selection of the type of data 
provided to the virtual identity system modeler 100 
5 and are not typically a physical electrical switch. 
Rather, switches lOlA and lOlB represent software 
routines or modules used to acquire and provide the 
selected data to the virtual identity system modeler 
100 . 

10 In FIG. 5, the inverse model 98 includes a 

static inverse model component 98A and an adjustable 
component 98B. The static component 98A is similar to 
the inverse model 72, discussed above. For instance, 
the static component 98A can be the inverse model 

15 H(f)."^ that was calculated by taking the inverse of 

forv;ard model H(f) at step 58 (FIG. 3A) . The static 
component 98A of the inverse model typically comprises 
a M X N matrix that includes cross-coupling effects, 
v/here M is the number' of inputs (drive 17) and N is 
^2 0" the numiDer^of' outputs' (actual re spons e 21 ) 7 ' 

The adjustable component 98B receives the 
inverse model correction values 104 from the virtual 
identity system modeler 100 for purposes of 
iteratively adapting the inverse model 98 to the 

2 5 current operating conditions. In one embodiment, the 

adjustable component 98B comprises an N x N matrix 
with . correction values for each of the N channels 
located on the diagonal and all other values (off- 
diagonal) equal to zero. 

3 0 FIG 6 is an embodiment of the functionality 

discussed in FIG. 5 in the context of the overall 
iterative control process detailed in FIG. 4B. 
Relative to FIG. 4B, in FIG 6, the inverse model 98 is 




SUBSTITUTE SHEET (RULE 26) 



r 



wo 99/38054 PCT/US99/0 1 233 

-17- 

extended to include the static component 98A and the 
adjustable component 98B of FIG. 5. 

FIG. 7 ill-ustrates a method of operation for 
the embodiment of FIG, G. At step 13 0, the inverse 
5 model 98 is initialized. This step includes 
initialization of the static component 98A as 
discussed above with respect to the method of FIG. 3A 
(step 58) , and initialization of the adjustable 
component 98B,. which for the first iteration is 
0 typically an identity model. In other words, the 
adjustable component 983 has no effect during the 
first iteration. 

At step 132, the initial drive 17 is 
obtained by convolving the target response 7 9 with the 
5 inverse model 98. In view that the inverse model 98 
includes the static component 98A and the adjustable 
component 98B, the target response 79 is first 
convolved with the adjustable component 983, the 
output of which is then convolved with the static 

1 comppnenjt;„_98A.. .In_certain. cases,, the static, component 

98A and the adjustable component 98B can be combined 
subsequently, requiring a single convolution. At step 
13.4, the drive 17 is applied to the physical system 10 
wherein the actual response 21 is measured and 
> recorded. 

Having obtained the actual response 21 of 
the physical system 10 from the initial drive 17, the 
virtual identity system modeler 100 performs spectral 
analysis between the target response 79 and the actual 
) -response 21. In the embodiment illustrated in FIG. 6, 
the virtual identity system modeler 100 includes a 
spectrum analyzer 13 6 that receives the target 
response 79 and the actual response 21. ' In one 
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embodiment on a channel by channel basis, the spectrum 
analyzer 136 calculates an FRF between the target 
.-response 79 and the actual response 21. This is 
represented at step 140. In other words, this step 
5 calculates a quality of the identity of the inverse 
model 98 and the physical system 10 (i.e. the virtual 
identity system 96) . The values 104 (FIG. 6) are 
corrections to the inverse model 98. At step 142, the 
deviation of the FRF from identity is determined and 
10 if one or more of the deviations exceed the 
corresponding selected threshold, values on the 
diagonal of the adjustable component 98B are updated 
at step 144. This is represented in FIG. 6 where 
-previous values of the adjustable component 98B stored 
15 at 14 8 are combined with new values provided by the 
spectrum analyzer 136, iteratively correcting the 
adjustable component 98B. 

At this point, it should be noted that steps 
88, 90, 97 and 99 of FIG. 2C are still performed in 
._2.Q. . order to obtain^ the. .refined- drive— 1-7 However, - before 
the new drive 17 is calculated at step 99, the inverse 
model is updated at step 144. 

In the overall iterative process, steps 134, 
140, 142 and 144 are only repeated as necessary when 

2 5 deviations in the model between the actual response 21 

and the target response 79 are greater than the 
selected model threshold. Steps 88, 90, 97 and 99 of 
FIG. 3C are performed independently for each iteration 
until the error threshold of step 90 is realized. 

3 0 It should be understood by those skilled in 

■ the art that although the virtual identity system 
modeler 100 and corresponding inverse model 98 are 
described in terms of spectral analysi-s methods, other 
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mathematical models and model regression techniques, 
either parametric or non-parametric, can be employed 
as desired in selected combinations (e.g. AIC, ARX, 
ARMA, State Space) . 
5 FIG. 8 illustrates an alternative embodiment 

where the actual response correction 103 is compared 
to the target response correction 77 and is used as a 
basis for updating the values of the adjustable 
component 98B. -In FIG. 8, a* summer 160 is used to 
0 obtain a difference between the actual response 21 and 
the immediately preceding actual response 87 (i.e. the 
actual response correction 103). During the iterative* 
process, the spectrum analyzer 136 calculates an FRF 
between the target response correction 77 and the 
5 actual response correction 103 in order to update the 
adjustable component 98B. 

In FIG. 8, target response correction 77 is 
added to the previous target response 79A to form the 
new target response 79 ..fsr each iteration. As 
, appreciated ,by those , skilled . in . the .art, _the target 
response correction 77 alone can be convolved with the 
.inverse model 9 8 to form the corresponding drive 
correction 94," which can be combined with the previous 
drive 17A to form new drive 17. Formation of the new 
drive 17 in this -manner does not alter the comparison 
by the spectrum analyzer 136 

Other general aspects of the present 
invention are illustrated in FIGS. 9 and 10. For cases 
where the target response 79 (or target response 
correction 77) does not exist in the forward iteration 
control loop, it can be explicitly computed with 
mathematical . equivalence as demonstrated in FIG 9 . In 
FIG. 9, drive 17 is applied to the physical system 10 
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and to a forward model 172 of the physical system 10, 
wherein a virtual identity system is indicated by 
-dashed lines 175.- It can be shown that virtual 
identity system 175 is mathematically equivalent to 
5 the virtual identity system 96 of FIG. 5. 



system 10 and a modeled target response 176 from the 
forward model 172 are provided to the virtual identity 
system modeler 100, In one embodiment, the virtual 
10 identity system modeler 100 performs spectrum analysis 
between the actual response 21 and the modeled target 
response 176 to check the quality of . identity of 
virtual identity system 175, since signal 176 is equal 
to either signal 79 or signal 77 of FIG. 5, depending 
15 on the position of switch 177A. The virtual identity 
system modeler 100 subsequently adjusts the model 172 
accordingly as a function of the quality of identity. 

In the embodiment illustrated, the model 172 
includes a static component 172A and an adjustable 
2-0-..- component 172B-. -The. .static- component 1-72A- can be- 
obtained per step 58 of FI.G. 3A. The static component 
172A typically comprises a N x M model that includes 
cross-coupling terms . 

The virtual identity system modeler 100 
25 provides correction values 104 to the adjustable 
component 172B. In one embodiment, the adjustable 
component 172B comprises a N x N diagonal model , 

In a manner similar to FIG. 5, switches i77A 
and 177B allow the inputs to the virtual identity 
3 0 system modeler 100 to be either the modeled target 
' response and the actual response or the modeled target 
response corrections and the actual response 
corrections. Those skilled in the art will recognize 



The actual response 21 from the physical 
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that switch 177A and associated summer and storage 
device could also be applied in the data path between 
the static model 1-72 A and the adjustable component 
172B rather than operating on the drive 17 as • 
5 illustrated. Likewise, the switch 177A and associated 
summer and storage device can also be applied to the 
modeled target response 176 in some circumstances. 

Alternatively, a similar implementation that 
could have preferred characteristics, potentially such 
10 as when N and M are unequal with the number of outputs 
or responses N greater than the number of inputs or 
drives M, is illustrated in FIG. 10. In FIG. 10, a 
virtual identity system is indicated by dashed lines 
185 and comprises physical system 10 and inverse model 
15 98 . . However, the quality of the identity is 
ascertained by comparing drive 17 applied to the 
physical system 10 with a corresponding modeled drive 
signal 182 obtained from the inverse model 98. Note in 
this case, the virtual identity system 185 is formed 
20. .with, respe^ct to. drive signals . as . opposed, to the use of 
response signals as in other embodiments. As always, 
the actual response 21 is obtained from the physical 
system 10 when the drive 17 is applied thereto. The 
actual response 21 is then provided as an input to the 
25 inverse model 98. In one embodiment, the virtual 
identity system modeler 100 performs spectrum analysis . 
between the drive 17 and the modeled drive 182 for 
successive iterations . The virtual identity • system 
modeler 100 subsequently adjusts the inverse model 98 
30 and, more particularly, the adjustable component 98B. 

In a manner similar to FIGS. 5 and 9, 
switches 183A and 183B allow the inputs to the virtual 
identity system modeler 100 to be either the drive 17 
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and the modeled drive 182 (as illustrated with 
switches 183A and 183B in the "G" position) , or the 
drive corrections 94- and the modeled drive corrections 
(as illustrated with the switches 183A and 183B in the 
5 "L" position) . Those skilled in the art will be able 
to recognize that switch 183B and associated summer 
and storage device can also be applied in the data 
path between the static inverse model 98A and the 
adjustable component 98B rather than operating on the 
0 actual response 21, as illustrated. Likewise, the 
• switch 183B and associated summer and storage device 
can be applied to the modeled drive signal 182 in some 
circumstances. 

FIG- 11 is an embodiment of the 
5 functionality discussed in FIG. 9 in the context of 
the iterative adaptive process of FIG. 8 where the 
same reference numerals have been used to identify 
identical components. FIG. 11 further includes 
relaxation gain 179 and frequency weighting function 
0., ,181 . Relaxation, gain 179 is. similar^ ,to„ relaxation, gain 
95, but provides gain on the drive correction 94 
rather than on the response error 89. Frequency 
weighting function 181 can be manually defined as well 
as computed from coherence type quantities (as 
indicated by "C" in 181) or other measures of model 
quality, generally as, but not limited to a function 
of frequency. An example of a coherence type quantity 
for frequency weighting of the drive correction- 94 at 
181 can be formulated as [H2] [HI] , where HI is a 

forward system model for the physical system 10 
assuming noise on the inputs, and H2 is a forward 
system model assuming noise on the outputs. Similarly, 
frequency weighting (e.g. [HI] * [H2] ""^ ) on the responses 
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can be provided at 18 9, if desired. 

The presence of either relaxation gain 179 
or frequency weighting 181, or other such functions, 
between the inverse model 98 and the phy3ical system 
5 10, breaks up the virtual identity system 96 of the 
forward iteration loop illustrated in FIG. 5. 
Therefore, in the embodiment of FIG. 11, the virtual 
identity system 175 is explicitly constructed per the 
illustration of FIG. 9. 
10 Referencing FIG. 11, the forward model 172 

includes the static forward model 172A and the 
adjustable component 172B of the physical system 10. 
Per step 58 of FIG. 3A, the static forward model 172A 
is determined; it's association with the physical 
15 system 10 is represented in FIG. 11 by dashed line 
190. With the static forward model 172A determined, 
the static inverse model 98A of the physical system 10 
is then calculated, the association represented by 
dashed line ^3 3oth the forward adjustable component 

2 0 172B-- -and - inverse-- - ad j us table - ~ component. 986^- are then 

initialized with diagonal elements set to one and off- 
diagonal elements set to zero. 

Initial drive 17 is obtained consistent with 
the method discussed above. The initial drive 17 is 
2 5 then applied to the physical system 10 wherein the 
actual response 21 is measured and recorded. As 
illustrated, the initial drive 17 is also applied to 
the forward model 172 to generate the modeled response 
176 (an equivalent response to the target response 79 
30 of FIGS. 6 and 8 with switches 177A and 177B in the 
"G" position) . 

The virtual identity system modeler 100 
compares the actual response 21 to the modeled 
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response 176. In the embodiment illustrated, the 
virtual identity system modeler 100 comprises the 
spectrum analyzer IS 6 and performs spectral analysis 
between the actual response 21 and the modeled 
5 response 176 on a channel by channel basis. Consistent 
with previous embodiments, the deviation of the FRF 
from identity is then determined by the spectrum 
analyzer 136 and if one or more of the deviations 
exceed the corresponding selected threshold, values on 

10 the diagonal of the adjustable component 172B are 
updated. By simply taking the inverse of the 
adjustable component 172B, values for the adjustable 
component 98B can be easily calculated for the next 
iteration, which begins with the calculation of the 

15 response error 89. 

A relaxation gain 187 can provide 
stabilization of the model update iteration loop. The 
rela:>f:ation gain 187 can be applied to correction 
values 104 as a -"power" k^, thereby creating relaxed 

'2'0*"^*'correctXo'fr values' 178 tKat" ar*e" appalled* Fo' ' the 'm 

correction 172B. Note, also, that the relaxation gain 
187 gain can be incorporated in any of the -embodiments 
discussed in this application. 

FIG. 12 is an exemplary embodiment for generating 

25 drive signals using the virtual identity system 
modeler 100 described above with respect to FIG. 10, 
where the same reference numerals have been used to 
identify identical components of previous embodiments. 
As in FIG, 10, the virtual identity system 185 is 

30 characterized using drive signals. The virtual 
identity system modeler 100 checks the quality of the 
identity by comparing drive 17 applied to the physical 
system 10 with a corresponding modeled drive 182 
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obtained from an inverse model 198. The inverse model 
198 is identical to the inverse model 98 and comprises 
the static component 98A and the adjustable component 
98B. The embodiment of FIG, 12 operates in a manner 
5 similar to the embodiments of FIGS. 6 and 8 wherein 
the virtual identity system modeler 100 (including the 
spectrum, analyzer 13 6) adjusts the inverse model 198 
by providing update values (herein relaxed correction 
values 178) to the adjustable component 98B thereof. 
0 However, in this embodiment, the adjustable component 
98B of the inverse model 98 is also updated in 
accordance with changes made to the adjustable 
.component 98B of the inverse model 198. 

FIG. 13 illustrates use of the virtual identity 
5 system modeler 100 in a Spectral Density Control 
embodiment. Unlike Time History Control that seeks to 
reproduce the response of the remote transducers 2 0 
with respect to time, Spectral Density Control seeks 
t'T reproduce the signal power in the response 

3 injc_luding„generat ing„ the__crossrpower .be.tween,, channels.) 

as a function of frequency over a selected bandv/idth. 
A power spectral density (PSD) comprises a square 
matrix with the auto power of each channel on the 
diagonal and the cross -powers between channels on the 
3 respective of f -diagonals . 

In FIG. 13 the same reference numerals have been 
used to identify similar components as described with 
previous embodiments. Generally, a summer 200 
calculates a PSD error 202 between a desired PSD 
) response 2 01 and an actual previous PSD response 203A 
from the physical system 10. The calculated PSD error 
202 is functionally similar to the time history 
response error 89 calculated in the previous Time 
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History Control embodiments. Commonly, a relaxation 
gain 204 is applied for substantially the same reasons 
as the relaxation gain' 95 to generate PSD response 
correction 206. In the embodiment illustrated, the PSD 
5 response correction 206 is combined with a previous 
PSD target response 207A to generate a new PSD target 
response 207. A PSD- to-time converter 208 converts the 
PSD target response to the time equivalent target 
. response 210 (sitnilar to target response 79) that is 
.0 provided to the inverse model 98. The inverse model 
98, in turn, is used to generate a drive 212 (similar 
to drive 17) that is applied to the physical system 
10, An actual response 214 (similar to actual response 
21) from the physical system 10 is provided to the 
5 virtual identity system modeler 100 (herein the 
spectrum analyzer 136) and to a time-to-PSD converter 
216. The virtual identity system modeler 100 checks 
the quality of identity in a manner similar to the 
embodiment of FIG. 6 and updates the inverse model 98, 
0_.. and ,in_ „ part icular the. ad j us table .component .98B 
accordingly. The time-to-PSD converter 216 generates 
the actual PSD response 203. As appreciated by those 
skilled in the art, any of the teachings of the 
previous embodiments of FIGS. 8-12 can be applied by 
5 replacing substructure 219 with the corresponding 
structure of .the previous embodiments. 

Another embodiment of Spectral Density 
Control is illustrated in FIG. 14, wherein the drive 
spectra are calculated directly from the response 
0 spectra in the forward loop as shown in block 18 9 (as 
familiar to those skilled in the art) , where the 
inverse system model is obtained from 172. The virtual 
identity system modeling concepts of FIG. 11 are added 
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to the forward spectral density iteration loop to 
iteratively adapt the system model 172 for the reasons 
discussed. 

FIG. 15 illustrates use of a Waveform 
5 Control embodiment. Like Time History Control, 
Waveform Control seeks to reproduce the response of 
the remote transducers 20 with respect to time. 
Waveform Control does so however without the feedback 
of the actual response into a summer to provide the 

10 response error 89. Rather, each iteration in Waveform 
Control uses the desired response 22 directly, applied 
through an attenuation factor 220 to make the target 
response 79 . 

Prior art techniques in Waveform Control 

15 recalculated the inverse model each iteration by using 
drive 17 and actual response 21 for that iteration in 
an effort to reach convergence. However, in this 
aspect of the present invention, the virtual identity 
zy^'cem 96 includes the inverse model 98 (static 

,2.0 component .98A„ and, adjus_table ,cpi^po^_^t: ?^8_B) a^nd 

physical system 10 as illustrated in FIG. 15. The 
virtual identity system modeler 100 receives the 
target response 79 and the actual response 21 to 
measure the quality of identity. As in the previous 

25 embodiments, the adjustable component 98B is updated 
as a function of the quality of identity. 
Incorporation of the virtual identity system 96 
improves single channel system operation and 
facilitates extension to multi-channel systems not 

30 previously feasible. 

FIG. 16 is a pictorial representation 230 of 
an example time history record of remote transducer 
response data from a statistically non.- stationary road 
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surface; where a first section 232 is indicative of a 
sequence of potholes, and where a second section 234 
-is indicative of cobblestones. Another aspect of the 
present invention includes constructing a virtual 
5 identity system for each section 232 and 234, and 
obtaining an adjustable component 98B for each of the 
sections 232 and 234 independently of a common static 
component 98A. Each adjustable component 9 8B is 
applied in combination v/ith the common static 
10 component 98A to obtain suitable drives for each of 
the sections 232 and 234, which are further combined 
to obtain drive 17 for the entire record 230. Note 
that virtual identity system modelers are not required 
to apply to contiguous sections of the record 230. 
15 Referring to FIG. 17, another aspect of the 

present invention includes allowing the model to vary 
within a section (e.g. section 234) or over the entire 
record 230, as illustrated below. Using this time 
varyipg^'-approach, the spectrum analyzer 136 operates 

20 on,, consecutive _ and., preferably., overlapping, analysis, 

windows 24 0 from each of the target and actual 
responses. Each of the target and actual responses is 
advanced by a selected time step 242, . thereby 
producing a sequence of spectral values 244 (as 
25 opposed to one net spectral average implied in the 
previous embodiments) . In this case, each of the 
target and actual spectral averages is formed as a 
two-sided running average of the respective individual 
spectral values, hence corresponding time varying 
30 spectral averages 246 result. The Step. 242 is 
typically between 10-90 percent of the analysis window 
or frame 240. Preferably, the step is 10-50 percent of 
the analysis window. 
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The individual elements of the time varying 
spectral averages 246 are processed to form a time 
varying model 104 (erg. FRF) and a time vairying model 
correction component similar to 172B per techniques 
5 discussed in the previous embodiments. Each correction 
component member 172B and corresponding inverse 
correction component member 98B of the sequence of 
correction components advances in respect to the time 
history, whereby processing occurs in time steps equal 
10 to time step 242. Generating drive 17 involves 
applying each correction component 172B and 
corresponding inverse correction component 9 8B of the 
sequence of the correction components to the 
respective time step 242 of the input record of .the 
15 target response 23 OB and combining the results thereof 
to produce a contiguous drive 17. 

Summarizing this embodiment, consecutive and 
preferably overlapping analysis windows 24 0 of the 
record 230 are formed in a stepwise fashion. In the 

--.2 0 context,- - of ^the spectral analyzer ^13 6,_. individual 

spectral values 244 are obtained for each analysis 
window 240. Individual spectral values 244 are the 
combined to form two-sided running spectral averages 
and corresponding FRFs , which are used in a step-wise 
25 manner to generate drive 17 or drive correction 94. 
Note that other statistical functions can be performed 
instead of or in addition to averaging. 

The spectral running average environment has 
an inherent stepwise character for processing the 
30 model update and applying the model correction. In an 
alternative embodiment, the correction model can be 
implemented in parametric form (i.e. not window based) 
with associated model regression method types such as 
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AIC, ARX, allowing the model to vary from sample point 
to sample point of the record 230. While the forward 
static model would be -physically realizable, the 
forward model correction is generally not physically 
5 realizable. Furthermore, because the adaptive process 
is occurring iteratively, at each sample point both 
past and future data are generally available to 
optimize the model regression process. 

FIG 18 is an embodiment of the present 
0 invention that facilitates adjusting the gain on the 
•time history input to the inverse model 98 on a sample 
point by sample point basis. As illustrated, sample 
point gain adjusting block 250 receives, as input, the 
target response 176 and the actual response 177B, or 
5 the .corrections thereof based on the position of 
switches 177A and 177B. Note these are the same 
fundamental inputs received by the virtual identity 
system modeler 100. 

The basic operation of the sample point gain 

.0 adjusting^blpck , 250. .is . to ^relate the., target .response., 

to the actual response (or corrections thereof) at 257 
such that a ratio or gain is realized on a sample 
point by sample point basis, indicated as a (k) , thus 
modeling the correction gain. The output gain per 
5 sample point is then applied at 254 to adjust the time 
history input to the inverse system model 98. 

In one embodiment, as indicated in FIG. 18, 
it may be desirable to apply filtering and threshold 
operations 256 and 258 to the input signals, which is 
0 illustrated as forming part of the sample point gain 
adjusting block 250. Similarly, it may be desirable to 
provide a filter 260 to filter the output of the 
sample point gain adjusting block 250.. As appreciated 
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by those skilled in the art, sample point gain 
adjusting block 250 and block 254 can be incorporated 
'in any of the Time History Control embodiments 
described above. 
5 Although the present invention has been 

described with reference to preferred embodiments, 
workers skilled in the art v/ill recognize that changes 
may be made in form and detail without departing from 
the spirit and scope of the invention. For example, in 

10 FIG. 11 (similarly in all embodiments) , for certain 
model types it can be advantageous to directly update 
(i.e. non-iteratively) the inverse model correction 
98B. In this scenario, the virtual identity system 
modeler 100 would receive as an input model response 

15 from static component 172A for all iterations. 
Further, the relaxed correction values 178 from the 
virtual identity system modeler 100 are directly 
applied as the inverse model correction 98B. 

-2 0- *. ■ - • • . 
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WHAT IS CIAIMED IS: 

1 . A method of controlling a physical system 
-responsive to an input to produce a selected output, 
the method comprising: 

defining a virtual identity system which includes 
the physical system for receiving the drive 
to provide an actual output; and 

checking a quality of identity of the virtual 
identity system using at least a function of 
the actual output . 

2 . The method of claim 1 wherein the virtual 
identity system includes at least one of an inverse 
model of the physical system and a forward ■ model, of 
the physical system. 

3 . The method of claim 2 wherein said at least one 
of the inverse model of the physical system and the 
forward model of the physical system includes an 
adjustable, component-, and- a. static^ component . . 

4. The method of claim 3 and further comprising: 
adjusting the adjustable component as a function 

of the quality of identity. 

5. The method of claim 2 and further comprising: 
adjusting said at least one of the inverse model 

of the physical system and the forward model 
of the physical system as a function of the 
quality of identity. 

6. A computer readable medium including instructions 
readable by a computer, which when implemented, cause 
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the computer to control a physical system responsive 
to an input to produce a selected output, the 
-instructions perfor^ning steps comprising: 

defining a virtual identity system which includes 
the physical system for receiving the drive 
to provide an actual output ; and 
checking a quality of identity of the virtual 
identity system using at least a function of 
the actual outout. 



7. The computer readable medium of claim 6 wherein 
the virtual identity system includes at least one of 
an inverse model of the physical system and a forward 
model of the physical system. 

8. The computer readable medium of claim 6 wherein 
said at least one of the inverse model of the phys.ical 
system and the forward model of the physical system 
includes an adjustable component and a static:- 
component. 

9. The computer readable medium of claim 8 and 
further instructions for performing a step comprising: 

adjusting the adj us table component as a function 
of the quality of identity. 

10 . The computer readable medium of claim 7 and 
further comprising instructions for performing a step 
comprising : 

adjusting said at least one of the inverse model 
of the physical system and the forward model 

of the physical system as a function of the 
quality of identity. 
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11. A computer implemented method of controlling a 
physical system responsive to an input to produce a 
selected output, the computer implemented method 
comprising : 

defining a virtual identity system which includes 
the physical system for receiving the drive 
to provide an actual output ; and 

checking a quality of identity of the virtual 
identity system using at least a function of 
the actual output . 

12 . The computer implemented method of claim 11 
-wherein the virtual identity system includes at least 
one of an inverse model of the physical system and a 
forward model of the physical system. 

13 . The computer implemented method of claim 12 
wherein said at least one of the inverse model of the 

phy s i c a 1„ . sy st em. . and^ t he„ forward, . mode l.^-of — the-., phy s i c al. 

system includes an adjustable component and a static 
component . 

14 . The computer implemented method of claim 13 and 
further comprising: 

adjusting the adjustable component as a function 
of the quality of identity. 

15. . The computer implemented method of claim 12 and 

further comprising : 

adjusting said at least one of the inverse model 
of the physical system and the forward model 
of the physical system as a- function of the 
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quality of identity. 



16. A system controller for controlling a physical 
system responsive to an input to produce a 
selected output, the system controller 
comprising : 

means for defining a virtual identity system 
which includes the physical system for 
receiving the drive to provide an actual 
output ; and 

means for checking a quality of identity of the 
virtual identity system using at least a 
function of the actual output. 

17. The system controller of claim 16 wherein .the 
virtual identity system includes at least one of an 
inverse model of the physical system and a forward 
model of the physical system. 

18. The system controller of claim 17 wherein said at 

■ - -least- one— of— t-he—i-nverse- -model- of- the- physica-l— system- - 
and the forward model of the physical system includes 
.an adjustable component and a static component. 

19. The system controller of claim 18 and further 
comprising: 

means for adjusting the adjustable component as a 
function of the quality of identity. 

20. The system controller of 'claim 17 and further 
comprising : 



means for adjusting said at least one of the 
inverse model of the physical system and the 
forward model of the physical system as a 
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f unction of the quality of identity. 
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